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Abstract 

This paper describes the extension of the Modified Observables technique in 
estimating simultaneously more than one Trilinear Gauge Couplings. The opti- 
mal properties, unbiasedness and consistent error estimation of this method are 
demonstrated by Monte Carlo experimentation using Ivjj four-fermion final state 
topologies. Emphasis is given in the determination of the expected sensitivities in 
estimating the A 7 — Ag\ and Afc 7 — Agf pair of couplings with data from the 183 
GeV LEPII run. 



1 Introduction 



It has been shown [|TJ, that by expanding the probability distribution function (p.d.f.) and 
keeping only linear terms with respect to the Trilinear Gauge Couplings (TGC's), one can 
build estimators (the Optimal Observables) which are linear functions of the couplings 
around the expansion point. Furthermore, this linear dependence can be easily evaluated 
by theory. An efficient estimation of the couplings can be performed by inverting these 
linear relations. Such an estimation has the same accuracy as the unbinned maximum 
likelihood technique. 

The method of the Optimal Observables has been extended |2| to incorporate the in- 
fluence of the detector effects to the measurement of the kinematical vectors. An iterative 
procedure has been also introduced to ensure the consistency and optimality of the tech- 
nique, independent of the choice of the parametric expansion point. In the same paper, 
the optimal properties of this (Modified Observables) method have been demonstrated 
for one coupling fits to the 172 GeV LEPII data. 

In the meanwhile larger data samples are available from the 183 GeV LEPII run and 
the application of the Modified Observable technique to a simultaneous estimation of two 
couplings is very relevant. 

This paper concentrates on the simultaneous estimation of two TGC's by employing 
phenomenological models || where two couplings could deviate freely from their Stan- 
dard Model (S.M.) values whilst certain constraints are imposed on the other couplings. 
This paper is dealing with WW events produced in e + e~ annihilation, where one of the 
Ws decays leptonically whilst the other decays in two jets. A large sample of 60000 
WW — > tvqq Monte Carlo (M.C.) events was used to evaluate cross sections and other 
statistics, as well as their dependence on the coupling values by the M.C. reweighting 
procedure [|J. These events have been produced either by PYTHIA [[J (employing only 
the CC03 production diagrams) or by EXCALIBUR || (full 4-fermion production) at dif- 
ferent coupling values and they have undergone full detector simulation by the DELSIM 
simulation programme. Moreover these events have been reconstructed and selected 
by the same analysis algorithms as the real data || || accumulated with the DELPHI 
|l(J detector. The background contamination has been simulated by the production of 



the physics channels || || (TT| which produce final state topologies indistinguishable from 
the signal WW events. 

This paper is organised as follows: the statistical technique and its asymptotic proper- 
ties are described in Section 2, whilst numerical results obtained by M.C. experimentation 
are presented in Section 3. Finally, Section 4 contains the comparison with other tech- 
niques and the conclusions. 



2 Modified Observables in Multi-Parametric Fits 

The present study is focusing on two parameter (TGC's) estimations but this analysis 
can be extended to any number of parameters in a straight forward way. 

The probability distribution function, with respect to the observed kinematical vector 
fl, is expressed 0] as a function of the two couplings a± and a 2 as 



1 



■/ CTtot (ai, a 2 J 
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where 

V is the true kinematical vector which describes the events 
e(V) is the efficiency of observing an event produced at V 

R(V; Q) is the resolution function, i.e. the probability the true kinematical vector V to 

— # 

be measured as f2 
da(v >i,Q2) j g differential cross-section 

dv 

Ctat(o!i,Oi2) is the total cross-section and 
Sij = J Cij(V)dV. 

In [[|, it has been shown that the Optimal Observables including detector effects, in 
the neighbourhood of the parametric point {a^a^}, are defined as the mean values of 
the following quantities: 

-j Sy 01 (V ] a°,a° 2 )e(V)R(V-,n)dV 

zuit: a,-, , qj 9 ) = - - — - — — 

fyoo(V;ala° 2 )e(V)R(V;tydV 

Zo (it, tt, , Oi = = — (Z) 

! yQ0 {V;ocla^)e{V)R{V^)dy 

where the functions y K ,\(V; a l5 a 2 ) are expressed in terms of the differential cross 
section coefficients as: 

y 00 (V- a°, a°) = c 00 (y) + c i(f )a° + c 10 (f)a° + c 20 (f)af + c 02 (f)af + c n (V)a%<4 
y 10 (V;a^a 2 ) = c 10 (V) + 2c 20 (V)a° 1 + Cll (V)a° 2 (3) 
yoi(^;a?,ag) = c i(V) + 2c 02 (V> 2 ' + c n (y)«? 

It has also been shown that the Optimal Observables are linear functions of the cou- 
plings a\ and a 2 in the neighbourhood of (a®, a 2 ) [] i.e.: 

J z k {VL- a?, a° 2 )P(Q; a u a 2 )dVt = J z k (Q; a?, a° 2 )P(Q- a°, a° 2 )dQ + 

J Zi(Q; al a° 2 )z k (n; a?, a° 2 )P(Q; a?, a° 2 )dQ - 
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a?, a!j)P(fi; a?, a!j)d« • 
* fe (fi; a?, a°)P(ft; a?, a°)eft2)] ■ (a, - a°) (4) 



lr rhis is easily proven by expanding (Fj) in a Taylor series around { a?, ajj} and evaluating the mean 
values of (0) ignoring higher than first order in a\ — a® and ai — a° terms. 
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where k=l,2 

Thus, given a set of N experimentally measured vectors Q n (n=l,. . . ,N) the left hand 
side of (H) can be approximated as: 

I z k ({l; a?, a°)P(Q; a u a 2 )d^l ~ - £ z fc (ft n ; a?, ar°) (5) 

n=l 

The right hand side of (f|) can be calculated using the theoretical expression of the cross 
section as a function of the couplings, provided that the resolution and efficiency func- 
tions can be parametrized analytically. Then, a simple inversion of the linear system of 
equations ([3j) results to an estimation of the coupling values with the same sensitivity as 
with the maximum likelihood technique. 

In practice, neither the efficiency nor the resolution function can be parametrized analyt- 
ically. However, it has been shown that a very succesfull approximative way of using the 
basic concepts of the Optimal Observables in one TGC parameter estimations [f2| exists. 
That is the Modified Observable technique, which in this paper is extended to more than 
one TGC simultaneous estimations. 

Following the same steps as in 0, the functional forms of Zi(ft; a?, a®) and 
^(fi; a?, 0!°) i n © are approximated as: 

/n o on 2/01,10 (A; aj, oP 2 ) 

These are very good approximations, as indicated in figure |] where the mean values 
of the quantities yoi ' 1 ? T ^'n 1 'n^ are compared with the quantities yoi ' 1 ^°'n 1 '^ , for several 
expansion points a and a®. These mean values have been evaluated by using M.C. events 
produced with coupling values a\ and ofy and being observed with kinematical vector Q 
corresponding to a bin of Wll,1 °l n 'n 1 'n^ ■ 

The functional form of 21,2(^5 «?, ot^) in (|6]) is independent of phase space and other 
multiplicative (e.g. Initial State Radiation) factors and it was calculated by using the 
ERATO []T3] four-fermion matrix element package by folding the kinematical information 
corresponding to the two hadronic jets. 

Instead of calculating the terms of the right hand side of (f|), the dependence of the 
mean values of Zi^iP-'t a i, 0*2) (i n the following called Modified Observables ) on the pro- 
duction values of the couplings has been evaluated by reweighted M.C. [|]] integration. 
Figure || shows the surfaces /i(ai, a 2 ; a , a°) an d /2(«i, «2! «?, a^) (i n the following called 
calibration surfaces), which express the dependence of the product of each Modified Ob- 
servable with the number of expected events for luminosity of 50.23 pb~ 1 , as a function 
of the coupling values, for three initial parametric points. These products (instead of the 
Modified Observables themselves) are going to be used as estimators of the couplings, 
gaining more efficiency by including the extra information of the total number of the 
observed events [HI. 



2 z\ 2(0; a?, otn) are defined in [01 as the mean values of the quantities ¥°hl2^J2i>22l corresponding to 

a yaa(V;a u 1 ,a'-J,) 

kinematical vectors V produced with coupling values aim a 2 an d being observed in the phase space 
element f2 • d£l. 
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The couplings are estimated by comparing the calibration surfaces to the experimental 
measurements, that is to the products of the measured values of the Modified Observables 
with the number of observed events, which are simply expressed as: 

N 

d 1 ( Q t?,a°)=532 1 (a;a?,aS) 

N 

«?>«!!) ( 7 ) 

1=1 

Such comparisons are shown in figures @ and (U) between a large independent set of M.C. 
events used as a data sample and three pairs of calibration surfaces (/i(ai, a 2 ; a?, a 2 ) and 
f 2(011, 0C2', a?> a< 2)) evaluated at three different expansion points {a®, a®}- In these figures 
the intersections of the calibration surfaces with the planes defined by the experimental 
measurements, di^a^a®) and G? 2 (a5,a 2 ), are also shown. It is worth noticing that the 
estimation, which is the common point of the pair of lines in figures ^p, |4|c and f|f, is 
independent from the expansion point. This fact reflects one of the basic properties of 
the technique to be globally unbiased. 

However, the evaluation of the estimation confidence intervals is more complicated, due 
to the statistical correlations between the calibration surfaces as well as between the mea- 
sured quantities di^a^a®) and ^(cti, a^)- 

The covariant matrices M(oei, 0:2; a i, 0^2) (expressing the statistical accuracy of the cali- 
bration surface evaluation at the expansion point { (a°, a°) } ) an d V(a1, a®) (which is 
the covariant matrix corresponding to the measured quantities di(a°, a 2 ) and d 2 (a?, « 2 )) 
are calculated from the kinematical vectors of the reweighted M.C. and real events re- 
spectively. 

Then, assuming gaussian errors, the probability that the selected event sample supports 
coupling values equal to a\ and cx.2, is given by the Likelihood function: 



1 1 -» -» T 

L = ^ • exp[-- (p(a>1, a° 2 ) - F(a h a 2 ; a?, a°)) • 

W~ l ■ oP 2 ) - F( ai , a 2 ; a?, a° 2 ))] (8) 

where the vector D(a\,a^), the vector calibration function F(ai, a 2 ; ot\, a°) an d the co- 
variant matrix W(ai, a 2 ; a°, a 2 ) are defined as follows: 

j\( n o fK _ / rfi(a?,a°) \ /on 

f(ai,a2 ' ai ' aj) - ^ /-,(<*!, a* a°,a°) J (10) 

W(ati, a 2 ; a?, a°) = M(a x , a 2 ; a?, a°) + V(a°, a°) (11) 

Maximization of with respect to «i and ct 2 , provides the estimation of the coupling 
values, whilst the confidence intervals are evaluated^ by the asymptotic gaussian properties 
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[logL max — 1.205) for 70% confidence intervals 
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of the estimation distribution |T3| . 

A set of M.C. events produced with Standard Model coupling values ( 6000 events at 
{A 7 = 0, Ag\ = 0}), was used as data sample to demonstrate the asymptotic properties 
of such estimations. The A 7 , Ag\ couplings were simultaneously estimated by maximizing 
the likelihood function of (|D and the estimated coupling values are shown as functions 
of the expansion point in figure (|5p. The fact that the estimations are close (within 
the statistical errors) to the true coupling values, for the whole region of the expansion 
points, emphasizes the optimal properties of the method. However, the optimal estimated 
error is achieved |]] at expansion points close to the estimated values, where the linear 
dependance of the Optimal Variables holds. This is shown in figure [|c where three 70% 
confidence limit contours corresponding to different expansion points are presented for 
comparison. Obviously the optimum estimated sensitivity is achieved in the case where 
a\ = d\ and a° = c?2> where d\ and &i are the estimated values. The fact that the above 
condition also guarantees a correct error estimation, is demonstrated in the next section 
by Monte Carlo experimentation. 



3 Numerical results 

A series of M.C. experiments were used to demonstrate the optimal properties of the 
Modified Observable technique when two TGC's are simultaneously estimated by fitting 
finite statistical samples. 

Fully reconstructed four fermion EXCALIBUR events, produced with S.M. coupling 
values, were mixed with background events to form data sets corresponding to the lu- 
minosity of 50.23 pb^ 1 accumulated by the DELPHI detector at ^/s ~ 183 GeV. Each 
set consisted of 82, 101 and 39 events in average with an electron, muon and tau lep- 
ton in the final state, respectively. The average background contribution to each of the 
above subsets were 8.0, 1.4 and 8.3 events. The specific event multiplicity of each data 
set was chosen to follow poissonian distributions. Another set of fully four fermion and 
background reconstructed events, produced and selected as described in Section 1, was 
used to calculate cross sections and probabilties as well as their dependence on the TGC's 
by reweighted Monte Carlo integration. In fitting the data sets, the (A 7 , Agf ) and the 
(Ak 7 , Agf) TGC schemes were used || and a simultaneous estimation of the free cou- 
plings was performed. 

In order to take into account the differences in the production dynamics, the selection 
efficiencies and the background contamination between the final states (ivjj, £ — fM,e, r) 
the measured vector -D(a°, a!]) was defined as follows: 

D(a°a°)-( £U4K« 2 °)-50 \. (12) 

Where i = 1, 2, 3 stands for the three lepton tags whilst B\ , B\ denotes the expected 
contribution of the background events to the measurement. 
Similarly the calibration surface vector was defined as: 



V 2^=1 J2\ a ^i a 2, Oil, a 2) 



(13) 
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The asymptotic property of the log likelihood ratio |TB| was used to demonstrate 



the unbiasedness of the proposed techniques. That is, the x 2 (n.d.f.=2) probability of 
obtaining the specific value of A, where 

T ( „,true „.true\ 

A = -2-log L( ^ ^\ ] (14) 
L{a u a 2 ) 

in fits of the data sets should follow an equiprobable distribution. 

Furthermore, the consistency in evaluating the error matrix of the estimated couplings 
{£) in every fit, is checked by using the other asymptotic property |13| of the likelihood 



estimations to be gaussian distributed around the true parameter values. Thus for an 
unbiased estimation of the central values and for a correct error matrix evaluation the 
quantity 5: 



.true 



5 = f ^ _ ^rue «2 " 4™ ) (15) 

should follow a x 2 (n.d.f.=2) distribution. This property is demonstrated by presenting 

the x 2 ( n -d.f.=2) probabilities to obtain specific 5 values in fitting the data sets. 

The above tests of A and 5 x 2 -probability distributions can be considered as extensions 

of the sampling and pull distribution tests respectively, commonly used in one parameter 

fits. 

Due to the limited number of the available M.C. events, only sixty independent data 
sets could be constructed. Although the number of the data sets is enough to show the 
optimal properties of the proposed technique, the bootstrap procedure Q [|14| has been 
also used to construct a large number of semicorrelated data sets. 

Results of estimating the (A 7 , Agf ) and (Afc 7 , Agf ) couplings with the Modified Ob- 
servable technique are shown in figure || In both TGC schemes the optimal properties 
of the technique in estimating central values and error matrices are obvious. Specif- 
ically the sixty completely uncorrelated samples produce x 2 ( n -d-f- = 2) probabilities 
distributed with mean values close to 0.5 and root mean squares close to l/\/l2, whilst 
the equiprobable behaviour of the x 2 { n -d-f- = 2) probability values obtained by fitting 
the bootstrapped samples is striking. 

The x 2 behaviour of the A and 5 quantities is further used to quantify the sensitivity 



of this technique. Indeed such property |L3j ensures that the estimated values {di,cf 2 } 



follow a two dimensional gaussian distribution with a covariant matrix which characterises 
the average sensitivity in estimating the couplings. The covariant matrix elements (i.e. 
the variances and correlations of the couplings estimations) are found by fitting a 2-dim 
gaussian to the estimated coupling values from the 60 independent sets. These average 
sensitivities are summarized in Tables 1 and 2 for (A 7 , Agf ) and (A&; 7 , Agf ) estimations. 
The same uncorrelated M.C. sets of events were treated as if they have been collected by a 
"perfect" detector and the two pairs of couplings were estimated by an unbinned extended 
likelihood fit as well as by the Modified Observable technique Q. The average sensitivities 

4 The bootstrap procedure advocates that one can select randomly Af events to form a set from 
a pool of /C available events, and repeat the random selection to construct many bootstrapped sets. 
The distribution of statistics evaluated from each of the bootstrapped set approximates well the true 
distribution, as long as K. is big enough compared to Af. 

5 The true kinematical vector V of each event was used to calculate the matrix element and the 
calibration surfaces. In the following, when V is used, the methods and their results will be named as 
"perfect" . 
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obtained from these estimations ("perfect" detector extended unbinned likelihood and 
"perfect" detector Modified Observables) are also shown for comparison in Tables 1 and 
2 where the equivalence of the Modified Observables to the likelihood fit is obvious. The 
loss of sensitivity in the case of a realistic detector is a natural consequence of the loss 
of information due to the imperfect experimental resolution. However, the consistent 
inclusion of the detector effects in the realistic case guaranties consistent central value 
and confidence interval estimation. It is also worth noticing that in the realistic case, the 
evaluated errors and correlations in every individual Modified Observable estimation are 
gaussian distributed with means very close to the average sensitivities, as it is shown in 
figure 0. 

As a last point, figure [8] shows the sampling, the pull and the error distributions of 
a single coupling (A/c 7 ) estimation. Similar estimations of the same coupling ||, using 
172 GeV data samples, have been found to exhibit non gaussian tails. However, it was 
advocated that with small event samples, where the statistical error is large compared 
to the linear part of the calibration curves, the evaluated error from the fits is expected 
to underestimate the sensitivity of the technique. Obviously such pathologies are absent 
when the relative statistical error is smaller, as in the case of the data sample accumulated 
at 183 GeV. 

4 Conclusion 

In this paper the Modified Observables technique || was generalized in order to be applied 
for a simultaneous estimation of two couplings by deploying the appropriate TGC scheme 
@. The technique, including the detector effects and the background contribution, was 
demonstrated to be asymptotically a consistent estimator. This consistency was also 
shown to be independent of the initial expansion values. However the optimal sensitivity 
is achieved for expansion points close to the estimated values of the couplings. 
The properties of the technique, when fitting finite size event samples, were investigated 
by M.C experimentation. Sets of M.C. events, of the same size as the data samples 
accumulated by each of the LEP experiments at y/s ~ 183 GeV, were fitted to estimate 
the {A 7 , A<jf} and {Ak 7 , Agf} pairs of couplings. The distributions of these estimations 
demostrated the optimal behaviour (unbiasedness, consistent error matrix evaluation) of 
the technique. Moreover a comparison with the unbinned extended likelihood results 
shows that the Modified Observable estimators are practically reaching the maximum 
sensitivity. These two methods are completely equivalent at the "perfect" detector case 
(tables 1 and 2). A deterioration of the sensitivity (up to 20%) when dealing with realistic 
detectors is due to the imperfect resolution of the measuring apparatus. 
A comparison || between the sensitivity of several multiparametric TGC estimators, 
which include detector effects, shows that the Modified Observables are equivalent to 
the Iterative Optimal Variables and Multidimensional Clustering techniques |l5j whilst 
outperform classical methods of one or two dimensional binned likelihood fits. 
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Method 


A 7 - Ag{ 








P 


"Perfect" Extended Likelihood 


0.21 ± 0.01 


0.20 ± 0.01 


-0.73 ± 0.06 


"Perfect" Modified Observables 


0.22 ± 0.01 


0.21 ± 0.01 


-0.74 ± 0.06 


Modified Observables 


0.25 ± 0.01 


0.23 ± 0.01 


-0.74 ± 0.06 



Table 1: Comparison of the statistical properties of the technique proposed in this paper 
with the unbinned extended likelihood estimations. 



Method 


A/c 7 - Ag{ 




CAfc 7 




P 


"Perfect" Extended Likelihood 


0.35 ± 0.03 


0.14 ± 0.01 


-0.22 ± 0.08 


"Perfect" Modified Observables 


0.38 ± 0.03 


0.13 ± 0.01 


-0.25 ± 0.09 


Modified Observables 


0.44 ± 0.03 


0.15 ± 0.01 


-0.28 ± 0.10 



Table 2: Comparison of the statistical properties of the techniques proposed in this paper 
with the unbinned extended likelihood estimations. 
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Figure 1: The mean values (see text) of the quantities ^ 10 |^"o'°oj [( a );( c )>( e )] an d 

[(b),(d),(f)l as functions of the ^#44^ and respectively. These 

mean values correspond to kinematical vectors produced with couplings: 
cni(= Agf) = 0, a 2 (= A 7 ) = in (a) and (b) 
«i(= A#f) = 0, a 2 (= A 7 ) = -1 in (c) and (d) 
aii(= AgJ 1 ) = 0, a 2 (= A 7 ) = +1 in (e) and (f) 
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Figure 2: The /i(cki, ol^, ofy, f 2(011, a 2] o^) calibration surfaces for several expan- 
sion points a?(= A^) = , a$(= A 7 ) = in (a) and (b), a?(= Ag%) = 0, a§(= A 7 ) = +1 
in (c) and (d), a?(= A^) = 0, a° 2 (= A 7 ) = -1 in (e) and (f). 
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Figure 3: The calibration surfaces fi(ea, a.2'1 a?, aQ), f 2(011, a 2', 0;°) as functions of 
the couplings cci = Agf, 012 = A 7 at the expansion point {a? = 0, a® = 0}. The 
horizontal shadowed planes correspond to the experimental measurements. The two lines 
representing the intersection of the calibration surfaces with the measured values are 
shown in (c). 
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Figure 4: The calibration surfaces fi(cxi, a 2 ; a®), f^i^ii a 2] «?, o^) as functions of the 
couplings cei = Ag^, a 2 = A 7 at the expansion points {a® = 0, a® = —1} in [(a),(b)] 
and {a® = 0, a® = +1} in [(d), (e)]. The horizontal shadowed planes correspond to the 
experimental measurements. The two lines representing the intersection of the calibration 
surfaces with the measured values are shown in (c) and (f ) for the two pairs of expansion 
points, respectively. 
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Figure 5: In [(a) and (b)] the estimated a\ = Agf and «2 = A 7 coupling values as 
functions of the expansion points are shown. In (c) the evaluated 70% CL. contours for 
a® = 0, a° = (solid points), a? = 0.2, a® = —0.6 (squares) and at{ = — 0.2, a® = 0.6 
(stars) are compared. 
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Figure 6: The distributions of x 2 (n.d.f.=2) probabilities in obtaining A [(a),(b),(e),(f)] 
and 5 [(c),(d),(g),(h)] values in {A 7 ,A«?f} [(a),(b),(c),(d)] and {A/c^Atf} [(e),(f),(g),(h)] 
estimations by the Modified Observables technique. The lines with slopes consistent with 
zero in (b),(d),(f) and (h) are first degree polynomial fits to the bootstrap results. 
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Figure 7: Confidence interval estimations with the Modified Observables technique. The 
distributions of errors [(a),(b),(d),(e)] and correlations [(c), (f)] in estimating the {A 7 , Agf} 
[(a),(b),(c)] and {Ak 7 , Agf} [(d),(e),(f)] pair of couplings. The data points correspond to 
the 60 independent data sets whilst the histograms to the bootstrap results. The arrows 
indicate the average sensitivities summarized in Tables 1 and Table 2. 
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Figure 8: The sampling, pull and error distribution of Ak 7 estimation with the Modified 
Observable technique. 
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